Observation of a resonant four-body interaction in cold cesium Rydberg atoms 
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Cold Rydberg atoms subject to long-range dipole-dipole interactions represent a particularly interesting sys- 
tem for exploring few-body interactions and probing the transition from 2-body physics to the many-body 
regime. In this work we report the direct observation of a resonant 4-body Rydberg interaction. We exploit the 
occurrence of an accidental quasi-coincidence of a 2-body and a 4-body resonant Stark-tuned Forster process in 
cesium to observe a resonant energy transfer requiring the simultaneous interaction of at least four neighboring 
atoms. These results are relevant for the implementation of quantum gates with Rydberg atoms and for further 
studies of many-body physics. 

PACS numbers: 32.80.Ee,34.50.Cx 
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The physics of atomic systems at low densities (n < 
10 13 cm -3 ) can generally be described in terms of the ac- 
tion of electromagnetic fields and binary (2-body) interac- 
tions. However, a number of interesting effects arise when 
few-body or many-body interactions come into play. Notable 
examples include: 3-body recombination QjJ-|3I] , leading to 
molecule formation in optical or magnetic traps; trimer pho- 
toassociation J3l; and Efimov physics, leading to trimers and 
more recently tetramer formation J^-d] . 

Cold, highly-excited (Rydberg) atoms |0- 11 ] are a promis- 
ing playground for many-body interactions, due to their strong 
and long-range interactions together with the long interaction 
times available in a cold sample. This was first revealed in 
studies 111 2i 1 1311 on the broadening of Rydberg energy trans- 
fer resonances. Subsequent work further explored the broad- 
ening mechanisms 11 1 lL Il44l8ll . the influence of the system 



dimensionality iflitl . and dephasing fl9ll . A renewed inter- 
est in cold Rydberg systems has recently been triggered by 
their possible use in quantum computation B20I1. Long-range 
van der Waals I2II1 or dipole-dipole I22I [23[| interactions al- 
low control of the atomic excitation in a given volume, the so- 
called blockade effect, enabling the implementation of quan- 
tum gates 11241 12511 . In this respect, a careful control over the 
number of contributing partners is necessary as 3-body and 4- 
body effects could be large enough to interfere with the com- 
putation process [26]. 

Another key feature of cold Rydberg atoms is the abil- 
ity to tune interactions by simply using an external electric 
field. In the process known as a Forster resonance, the en- 
ergy of the final many-body state can be Stark tuned into 
resonance with the initial state, leading to a resonant energy 
transfer (FRET) H271 12811 . In this case, modeling the system 
requires including multiple atoms and solving the full three- 
dimensional many-body wave equation lfl8il29ll . However, the 
number of atoms which must be included for accurate results 
is still unclear B30U32I1 . It is therefore important to study cases 
where a small number of atoms play a dominant role. This is 
the case in the striking experiments on 3- or 4-body recombi- 
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FIG. 1. (Color online) Calculated four- atom energy for the two rel- 
evant 2-body resonances (ai) and (a 2 ) and the p — > d 4-body reso- 
nance (b) as a function of the electric field. 



nation where the few-body process is characterized by 

a resonant loss mechanism. 

In this Letter we present first results on a 4-body energy 
transfer process due to a 4-body Stark-tuned FRET resonance 
occurring in cesium at an electric field of 79.99 V/cm: 



4 x 23p 3/2 -> 2 x 23s + 23p 1/2 + 23d- 



5/2- 



(1) 



We take advantage of two nearly resonant Cs 2-body Stark- 
tuned FRET resonances, 



23p 3/2 + 23p 3/2 
24s + 24s 



23s + 24s (2) 
23j5 1/2 + 23d 5/2 , (3) 



occurring respectively at 79.94 V/cm and 80.42 V/cm. All 
states have | to ,7 | = 1/2 unless specified otherwise. For conve- 
nience, the states 23p 3 / 2 , 23s, 24s, 23pi/ 2 and 23ds/ 2 are la- 
beled as p, s', s, p', and d, respectively. The three resonances 
are illustrated in Fig. [T] and are denoted as p — > d (Eq. (Q~|)), 
p — > s (Eq. (|2]l) and s — > d (Eq. (O). Observing d popula- 
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FIG. 2. (Color online) Thirty-shot-averaged normalized time-of- 
flight signals for laser excitation to the Stark-mixed d, s, p' and p 
states. The four traces have been vertically offset for clarity, with the 
respective zero levels shown as horizontal dashed lines. The electric 
field is set to 79 V/cm, far from any FRET resonance. The gates for 
the d, s, and p signals are shown as gray shaded regions. 

tion after exciting the p state is the signature that four atoms 
exchanged energy. 

The experimental setup consists of a standard Cs magneto- 
optical trap (MOT) at the center of four parallel 60 mm by 
130 mm wire mesh grids of 80 /im thickness and 1 mm grid 
spacing. The center pair of grids is spaced by 1.845±0.01cm, 
and the outer grids are 1.5 cm from the inner grids. Voltages 
up to ±5kV can be applied arbitrarily to the four grids. Six 
additional small electrodes surround the excitation region at 
the grid edges to cancel stray fields. The central grid spacing 
has been calibrated by measuring Stark shifts of the 23p and 
22d states for fields up to 100 V/cm and comparing them to 
theory §3ll . 

The trapped atoms are excited to a Stark-shifted Rydberg 
state, nl, using the MOT light and two additional lasers via 
the scheme: 6s — > Qp — > 7s — > nl. The 6p — > 7s step uses 
a 1470 nm, lOmW laser frequency-locked on a Doppler-free 
feature in a Cs cell excited by resonant 6s — > Qp light. In order 
to avoid perturbations of the atomic cloud, the 1470 nm laser 
is switched on for only 500 ns, at a 10 Hz repetition rate, by an 
acousto-optical modulator (AOM). The first order AOM out- 
put is focused to a 400 /im diameter spot at the atomic sample 
and its intensity is chosen to avoid power-broadening the tran- 
sition. A cw Tksapphire ring laser, providing roughly 1.8 W 
at rs 795 nm, drives the 7s — > nl transition. This laser is 
switched in the same manner and the beam is focused to a 
300 /im spot diameter and perpendicularly overlapped with 
the 1470 nm beam in the atomic sample. The simultaneous 
500 ns pulses are short enough to avoid excitation blockade. 
We thus excite a Gaussian cloud of up to 2 x 10 5 p atoms, 
260 /im in diameter, with a maximum peak density, p p , of 
9xl0 9 cm" 3 . 

Selective field ionization is used to measure the populations 
of the various Rydberg states. A voltage ramp is applied to one 
of the grids, 1.5 /is after the beginning of the laser excitation, 
rising to 4.3 kV in 4/xs. The resultant ions are detected by a 
microchannel plate (MCP) detector, 210 mm from the center 



of the trapped cloud. The amplitude of the field ionization 
pulse is chosen to optimally isolate the d time-of-flight (TOF) 
from the other signals, as displayed in Fig. 12 where each state 
has been excited directly in a non-resonant electric field. We 
use the time gates shown to measure the population of each 
state. 

It is important in this experiment to retrieve the correct pop- 
ulations of the p, s, and d states from a single TOF signal. 
As one can see in Fig. |2] we have the ability to differenti- 
ate these Rydberg states, but it is also apparent that they are 
not cleanly separated into their respective gates, for example 
due to blackbody radiation or multiple ionization pathways. 
We therefore use traces, shown in Fig. 12 where each state 
is separately excited and the electric field is not resonant, to 
quantify the crosstalk between the various gates. We also ac- 
count for the detection efficiency of each state, using a volt- 
age which ionizes completely, and separately measuring the 
global MCP detection efficiency. We notice a slight overesti- 
mation of the measured cross-talk coefficients leading to neg- 
ative s or d population for small p excitation. We correct these 
coefficients and these corrections are included in the error bars 
shown in the experimental figures below. Throughout the pa- 
per, the recorded signals are converted to the actual number of 
atoms in each state via the matrix, 

d\ ( 2.016 -0.0644 -0.082\ (d\ 
s = -0.100 4.645 -0.275 s , (4) 
p) V 0-083 -3.147 4.149 / \pj gatc 

whose off-diagonal elements characterize the crosstalk. For- 
tunately, the d-state signal has minimal crosstalk with the oth- 
ers, and a more advanced signal analysis [ 34] is not required. 

Although 2-body Stark-tuned FRET resonances have been 
observed previously Il2[ 12211 . we present them below for com- 
pleteness. We first verify the 2-body p — > s resonance, as 
shown in Fig. [3^. The resonance is observed at 79.94 V/cm, 
and the flat-top resonance shape suggests that the transition 
is saturated. We observe a small field inhomogeneity in 
our experimental region, which we estimate to be around 5 
V/cm/cm, via broadening of the p laser excitation line. This 
corresponds to ±0.05 V/cm in the excitation volume, which is 
thus the electric field measurement resolution. Stark field laser 
excitation of the s state lets us measure the 2-body s — > d res- 
onance, as shown in Fig.[3t>. We can excite up to 1.4 x 10 5 24s 
atoms at a peak density, p s , of 8 x 10 9 cm~ 3 . The s — >■ d reso- 
nance is observed at 80.42 V/cm. We have used the Cs energy 
levels published by Sansonetti [35] for our calculations. Un- 
fortunately, the 24p 3 / 2 state energy uncertainty (0.03 cm -1 ) 
creates a large uncertainty in the position of the 2-body reso- 
nances; ±0.1 V/cm for the p — > s transition, and ±1 V/cm for 
the s — > d transition. While the positions of the resonances 
are not well defined, the slopes shown in Fig. Q] are stable to 
within 0.1%, and we have therefore aligned the calculated 2- 
body resonances with our experimentally measured positions. 

Once the 2-body resonances have been measured, the loca- 
tion of the 4-body resonance is accurately known. We observe 
greater than 1% p — » d transfer with a peak at 79.99 V/cm, as 
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FIG. 3. (Color online) Number of detected Rydberg atoms as a func- 
tion of applied electric field, a) percentage of detected s atoms when 
exciting the p state, b) percentage of d detected when exciting the s 
state, c) percentage of d detected when exciting the p state. The error 
bars account for the corrections to the cross-talk coefficients, which 
are larger for s than for d, and for the observed field inhomogeneity 
that affects more the p excitation, leading to larger error bars than in 
the s excitation case. The results of our calculation for four equidis- 
tant atoms are overlaid as dashed blue lines, with the p — > s calcula- 
tion amplified by a factor of two and the s — > d divided by a factor of 
three to coincide with experimental results. The three resonant field 
values are illustrated by the vertical dotted lines. 
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FIG. 4. (Color online) Figure a) shows the number of detected d 
atoms as a function of the applied electric field for exciting the p (red 
+) and s (blue x) states, with an s density comparable or greater 
in the s excitation than in the p excitation. The on-resonant 4-body 
process creates more than a factor of four more d atoms than the off- 
resonant s — > d 2-body process at 79.99 V/cm. Figure b) shows the 
number of detected d atoms when exciting the p — > s resonance for 
the \mj\ = 1/2 (79.94V/cm) and \mj\ = 3/2 (88.14 V/cm). We 
observe no significant d state population in the \mj\ = 3/2 case. 



shown in Fig. [3};. While the 4-body resonance partially over- 
laps the p h s 2-body resonance in field, it is important to 
recall that the d state signal is well separated in the TOF. 

To provide insight into the 4-body resonance, we have de- 
veloped a minimal toy model with four equidistant atoms ar- 
ranged as a tetrahedron. The four possible states \pppp) (ini- 



tial state), \ss'pp), \ss'ss'), and \ds'p's') (detected state) are 
coupled by dipole-dipole interactions, calculated between the 
in-field eigenstates of the Rydberg atoms 13311 . The final pop- 
ulations, shown in Fig.[3]as the blue dashed curves, are calcu- 
lated using the density matrix and the experimental peak den- 
sity and field inhomogeneity. We average the results assum- 
ing an Erlang (nearest neighbor) distribution for the 2-body 
distance between the atoms, and a cubic Erlang distribution 
for the 4-body case. Such a 2-body model is not expected 
to precisely match the experiment [18] and the calculated 2- 
body p —> s curve is amplified by a factor of two while the 
s — > d curve has been diminished by a factor of three to match 
the experimentally observed results. To account for the 0.3% 
background observed in the p — s> d 4-body transfer, the 4-body 
curve baseline has been shifted accordingly. While a more de- 
tailed many-body calculation would be needed to reproduce 
the data 1I29I1 . it is remarkable that such a crude 4-body calcu- 
lation qualitatively reproduces the shape of the experimental 
signal. 

The observation of d state population constitutes a clear, 
direct signature of an interaction involving at least four bod- 
ies. The strong signature that the process is not a simple com- 
bination of two consecutive 2-body processes, but a genuine 
4-body process, lies in the relative strengths of the d transfer 
at the 4-body resonance field when initially exciting the p or 
s state. Figure |4^ shows the number of detected d atoms for 
comparable densities of s, either excited directly (1.4 x 10 5 ) 
or obtained from exciting p and allowing 2-body transfer into 
s (< 10 5 ). The larger number of detected d atoms (about a 
factor of four) when exciting p, despite a smaller s density, 
is explained by the fact that here the 4-body p — > d transfer 
is resonant whereas the 2-body s — > d transfer is not. Fi- 
nally, exciting the p state in the \mj\ = 3/2 case, where the 
p — > s resonance lies around 88.1 V/cm and the 4-body reso- 
nance is well separated from both 2-body resonances, we see 
in Fig. |4j5 that no significant d population is detected at the 
p H> s resonance. The small observed signal is compatible 
with the estimated error on the inversion matrix coefficients 
in Eq. (@). 

A way to increase the population transfer is to shift the ap- 
plied field from the p — > s resonance to the s — > d reso- 
nance between laser excitation and detection. The two 2-body 
FRETs are then consecutively resonant and we should get a d 
transfer comparable to that obtained when exciting s directly. 
Indeed, we have observed up to 7.5% p —> d population trans- 
fer with a 0.6 V/cm field shift, confirming that the p — > s pop- 
ulation transfer at 79.94 V/cm leads to about the same s den- 
sity as direct s excitation. With a larger shift it is also possible 
to induce a transfer to d from the \mj\ = 3/2 p state. We ob- 
served up to 2. 1 % p — > d population transfer with a -7.7 V/cm 
field shift starting from 88.14 V/cm. 

The field inhomogeneity and the proximity of the 4-body 
and 2-body FRET resonances impede further studies on the 
resonance shape. Nevertheless, having identified the 4-body 
resonance, we can study the transfer dependence on the initial 
p Rydberg atom density. We vary this density by attenuating 
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FIG. 5. (Color online) Number of detected Rydberg atoms as the 
intensity of the Rydberg excitation laser is varied, a) Number of 
d atoms as a function of s detected atoms, b) Number of detected s 
atoms as a function of the total number of detected Rydberg atoms, c) 
Number of d atoms created vs the total number of detected Rydberg 
atoms. Best fits below 1.6 x 10 atoms (demarcated with a vertical 
dotted line) are shown as the dashed blue lines. 



the Rydberg excitation laser with a set of neutral density fil- 
ters, sitting on the 79.99 V/cmp — > d resonance. The results 
are shown in Fig. [5] We have fit data only above the detection 
sensitivity of 500 atoms for s and 100 atoms for d. The popu- 
lation transfer from the initial p state to the s state is shown in 
Fig. [5J5. It increases as p 2 v until the p — > s transfer starts sat- 
urating above 1.6 x 10 5 atoms, as expected from Fig. [3] and 
we have limited all three fits to below this value. The d state 
population transfer as a function of the detected s population 
is shown in Fig.[5k. The quadratic, i.e. nonlinear, dependence 
ensures us of the correct data treatment presented in Eq. 
Finally, we observe in Fig. [5}: the number of detected d atoms 
as a function of the total number of excited Rydberg atoms. 
We clearly see the influence of the s transfer saturation on the 
d transfer. Nevertheless, we can see the d population trans- 
fer scaling as p*, again demonstrating that d state observation 
links to a 4-body process. 

In conclusion, we have observed a 4-body interaction be- 
tween Rydberg atoms excited in a MOT. The observed pro- 
cess is a FRET involving the simultaneous interaction of at 
least four neighboring partners. Its occurrence depends on the 
electric field proximity of the parent 2-body FRET resonances 
and to the quasi-coincidence of the 4-body resonance with one 
of them. We have studied the 4-body reaction as a function of 
the electric field and the initial Rydberg density by selectively 
detecting the reaction products. 

Tunable FRET resonances are extremely interesting tools 
for making Rydberg atoms into suitable systems for quantum 
information storage and processing 111 ll 12011 . Nevertheless, 



care should be taken to properly evaluate the occurrence of 
higher-order (3-body, 4-body or n-body) processes when is- 
sues of fidelity and decoherence are to be addressed. 

This work on a strong 4-body interaction is an important 
step towards studying few-body and many-body effects in di- 
lute gases. Such processes which go beyond the 2-body inter- 
action are likely to occur for many Rydberg states, especially 
at the larger densities that can be obtained in optical traps. For 
instance, other 4-body FRET resonances exist for four non- 
identical Rydberg atoms, for which the resonance field is bet- 
ter separated from the 2-body resonance field. Furthermore, 
the use of optical lattices to induce spatial order in the system 
should allow new insights into the physics of novel quantum 
systems. 
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